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SCALABLE DiMENSIONLESS ARi=lAY 
TECHNICAL FIELD 

This Invention relates to the general field of digital computing and in particular 
to a scalable array of globally docked multiply/accumulate floating point 
5 processing elements. 

BACKGROUND ART 

Kung and Leiserson ['Systolic Arrays (for VLSI)' in Sparse Matrix Proceeding 
1978, Soc. for Industrial and Applied Mathematics, 1979] presented the concept 
of performing matrix operations using arrays of simple processing elements. 
1 0 Each processing element implements a simple primitive operation. As an 
example, at a given time, a processor may: 

read the input data vector {a(in).b(in),c(in)}, 

perform an arithmetic operation such as c(out) = a(in)b(in) + c(in), 

write the output data vector {a(out),b(out),c(out)). 

15 The processing elements are connected only to their nearest neighbours, and 
so the problems of routing, fan-out and clock skew are minimised. Data and 
results move synchronously ttirough the array of elements. The name applied 
to this approach to computation with arrays of identical processing elements is 
systolic. 

20 An example algorithm quoted by Kung and Leiserson was the matrix product 
Using a systolic anay in which the processing elements executed the local 
algorithm presented above, known as an inner-product-step algorithm, they 
showed that the system level algorithm which this implemented was a matrix 
product of computational order 0(N), rather than the computational order 0(N3) 

25 for the matrix product implemented on a conventional scalar architecture. The 
matrix product is represented simply as 

C = AB 

where A, B and C are matrices of a size equal to the order of the anay. 
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Conformal matrices can be multiplied with ttiis array under certain restrictions if 
tlie results are re-drculated througti the array, and larger order matrix products 
can be computed if the tasic is partitioned. Speiser, Whitehouse and Bromley 
['Signal Processing Applications for Systolic An^ays', Record of the 14th 

5 Asiiomar Conference on Circuits, Systems and Computers, IEEE No. 
80CI-I1625-3, 1980 ] subsequently demonstrated the use of inner-product- 
accumulate processing elements for the same matrix product algorithm. In this 
case, the results are formed in-place, and do not move between processing 
elements. The only difference between the description of this algorithm and the 

1 0 algorithm described above is that the input and output phases of the algorithm 
do not include the reading and writing of c(in) and c{out) respectively, and that 
an explicit unload phase must be added at the end of the algorithm to return the 
results. 

The primary advantage of systolic processing over conventional linear 
1 5 processing is speed. The systolic architecture uses the fact that for matrix 
multiplication, the same operand datei may be reused many times in the 
computation of cross-product temis, thereby making better use of the available 
data bandwith. The Improved performance, however, comes at the cost of 
flexibility. Prior art devices have been designed for very specific applications 
20 such as Fast Fourier Transform computations or video signal processing. An 
advantage of the present invention is the ability of the same device to be useful 
for a wide variety of matrix computations without the need for hardware 
reconfiguration. The device is particulariy useful when implemented as an 
architectural enhancement to a computer in which case the processing power 
25 of the computer is considerably enhanced. 



DISCLOSURE OF THE INVENTION 



it is an object of this invention to provide a processing element for use in a 
scalable array processor which is able to Implement a set of primitive floating 
point matrix operations for conformable matrices of arbitrary order on an array J 
30 of fixed size. , 

It is a further object to provide a scalable array processor chip which is able to 
perform one or more of ttie following functions : 
compute the product of two matrices 

compute the element-wise (Hadamard) product of two matrices 
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compute the sum of two matrices 
pemiute the rows and columns of a matrix 
transpose a matrix. 

it Is a still further object of this invention to at least provide the public with a 
5 useful alternative to existing systolic devices. 

Therefore, according to perhaps one form of this invention, although this need 
not be the only or indeed the broadest form, there is proposed a processing 
element suitable for use in a scalable anay processor comprising : 
at least one input register means adapted to receive and process serial 
1 0 operands in the form of {instruction, data} 2-tuples; 

a memory means adapted to store temporary results and constants; 

a computing means adapted to perfomi logical operations; 

an output register means adapted to output results from the processing 

element; 

15 a control and sequencing means adapted to control the operation of the 
processing element; 

a plurality of data buses adapted to provide communication between the 
plurality of means. 

In preference the computing logical means consists of a shifter/normaliser 
20 means adapted to shift/normalise data and an arithmetic means adapted to 
perform logical operations such as but not limited to addition, subtraction and 
partial multiplication operations. 

In preference the processing element is adapted to perform floating point 
multiply, floating point add and floating point muttiply-accumulate which is used 
25 for inner product accumulate operations. 

In preference the input register is adapted to output a copy of ttie input operand 
bit with a one clock period delay. 

In preference tiiere are N input registers and the processing element Is suitable 
for use in a N-dimensionai scalable array processor. 

30 In preference N can be any positive integer. 
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in preference the input registers convert the Input serial data to an internal 
representation comprising separate sign, fraction and exponent 

In preference the memory means consists of read only memory for storage of 
constants and a read/write memory for storage of temporary results. 

5 In preference the shifter/nonmalizer means is adapted to perfonn binary 
weighted barrel shifting wherein the shifter function is detemriined by a control 
Input to the shifter/normallzer and the normalizer function effects a data 
dependent shift of up to 15 bits within a single clock cycle. 

In preference the arithmetic means implements logical operations such as but 
1 0 not limited to floating point addition, multiplication and multiply-accumulate 
algorithms using a parallel microcoded data path. 

In preference the arithmetid means comprises a logical unit such as but not 
limited to an input-multipiexer, an adder, an output shifter, flags unit and a 
control unit. 

15 in preference the output register can be loaded in parts to enable the 

conversion from the internal representation to IEEE 754 floating point format. 
The output register can be parallel loaded from the arithmetic means or can be 
serially loaded from a serial source. The register is unloaded serially. 

in preference the control and sequencing means includes timing and control 
20 logic, a microcode ROM, address decoders, branch control logic, flags logic, 
instoiction register, Instaiction decoder and a program counter. 

In preference there are three data buses, an X bus, a Y bus and a R bus. The X 
and Y buses are called operand buses and the R bus is called the result bus. 

In preference the processing element has an accumulator comparison means. 



25 In another form the invention consists of a scalable array processor chip 
comprising an array of processing elements each said element including : 
at least one input register means adapted to receive and process serial 
operands in the form of {instruction, data} 2-tuples; 
a memory means adapted to store temporary results and constants; 
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a shifter/normalizer means adapted to shift or normalize data; 
an aiitlimetic means adapted to perfomi logicai operations such as but not 
limited to addition, subtraction and partial mulfipiication operations; 
an output register means adapted to output results from the processing 
5 element; 

a control and sequencing means adapted to control the operation of the 
processing element; 

a plurality of data buses adapted to provide communication between the 
plurality of means; and 
1 0 wherein each processing element has means for communication only with 
adjacent elements. 

In preference the array of elements comprise an interconnected lattice of at 
least one dimension. 

« 

In preference the anay of elements comprise an Interconnected lattice of at 
1 5 least two dimensions. 

In preference the scalable anay processor chip is adapted to perfomi at least 
the functions of computing the product of two or more matrices, computing the 
element-wise product of two or more matrices, computing the sum of two or 
more matrices, permuting the rows and columns of a matrix and transposing a 
20 matrix. 

In a yet further form of the invention there is proposed a computing apparatus 
comprising a host processor, at least one scalable array processor chip and a 
plurality of data formatters wherein the scalable anay processor chlp(s) and 
plurality of data fomiatters are adapted to perfomi matrix operations othenvise 
25 performed by the host processor. 

In preference the apparatus includes a memory cache adapted to store 
operand data and temporary or intermediate resulte. 

In a still further form of the invention there is proposed a method of performing 
matrix operations comprising the steps of : 
30 (a) providing a plurality of processing elements in the forni of an array 
adapted to perform systolic processing operations; 

(b) receiving operand matrix data for processing from a host or data source; 
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(c) formatting the operand matrix data in a data formatter by adding an 
instruction to form an {instruction, data} 2-tuple; 

(d) transferring sets of 2->tuples to tlie processing element array to cause tiie , 
processing elements to process the data in accordance with the instruction; 

5 (e) repeating the steps (b) to (d) a number of times said number of times 
being dictated by the matrix operation being perfomied; 
(f) unloading the results of the matrix operations into an output result 
register of the processing elements (under control of an instruction specified by 
the operand 2-tuple); 
1 0 (g) transferring the contents of the output result registers held within the 
plurality of processing elements back to data formatters as result wavefronts; 

(h) storing the result wavefront data back to a host or data sink; and 

(i) repeating the steps (f) to (h) a number of times said number of times 
being dictated by tiie matrix operation being performed. 

15 in preference, the data formatter is of the type the subject of co-pending patent 
application number PL5696 entitied "DATA FORMATTER". 

The sets of 2-tuples are known as operand wavefronts. During the unload step, 
an indication is provided to the data formatter that the unload operation is 
occurring, allowing synchronization of data ti^ansfers to and from the processing 
20 array. 

During step (g), ttie results are transmitted in sets containing one resutt from 
each of the processing elements at the left edge of the array. Such a set is 
known as a result wavefront. 

There may be as many resutt wavefronte held wittiin the array as tiiere are 
25 columns of processing elements. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of this invention a preferred embodiment will now be 
described witii reference to the attached drawings in which : 

FIG. 1 is a schematic diagram of one processing element; 

30 FIG. 2 is a schematic diagram of one embodiment of a systolic array 
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processing element chip utilising the elements of FIG. 1; 

FIG. 3 isaschematicdiagramof a first embodiment of a processing 
apparatus utilizing the chip of FIG. 2; 

FIG. 4 is a schematic diagram of an inner-product-step processor; 

5 FIG. 5 is a schematic diagram of an inner-product-accumulate processor; 

FIG. 6 is a schematic example of the entry of operand wavefronts to a 
processor array; 

FIG. 7 Is a schematic example of the unloading of result wavefronts from 
a processor array; 

1 0 FIG. 8 Is a schematic example of tiie entry of element-wise operand 

wavefronts to a processor array; 

FIG. 9 is a schematic example of tiie unloading of element-wise result 
wavefronts from a processor array; and 

FIG. 10m is a schematic diagram of a second embodiment of a processing 
15 apparatus. 

BEST MODE FOR CARRYING OUT THE INVENTION 

Referring now to tiie drawings in detail, each processing element consists of a 
number of input registers, a memory consisting of a register file and a constant 
ROM, a shifter/nomializer, an arithmetic unit, output registers and a control and 
20 sequencing unit. 

The datapath elements ( input registers, memory, shifter/normaiizer, arithmetic 
unit and output registers) are interconnected by three parallel data buses. In 
addition serial interfaces are provided to and from each of the input registers 
and the output register to allow communication between processing elements 
25 and to facilitate construction of arbitrarily large arrays of processing elements. 

An array computes 2N2 floating point operations (1 multiply and 1 accumulate 
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for each processing element in the array) in the time taken to fetch 2N 
operands. To ensure that the computation is bandwidth limited, each 
processing element needs only compute at a rate of one floating point operation 
every N data fetches. This fact leads to the conclusion that very cheap 
5 processing elements can be used in the an^ay. 

A schematic of the processing element is shown in FIG. 1 . The choice of a 
simple microcoded datapath and sequential algorithms to perform the floating 
point operations means that the size of the processing element can be kept 
small. Many such processing elements can therefore be placed on a single 
10 chip. The fact that processing elements implemented in this manner are slower 
than those built using fully parallel algorithms and architectures becomes, 
insignificant as the size of the anray is Increased. This is because the 
processing performance achieved is limited by the data bandwidth (and.array 
size)t not by the computation rate for a single processing element 

1 5 The functions performed by each module in the processing element are 
described below. 

Input Registers: The input registers receive serial operands in the form of 
{instruction, data} 2-tuples from adjacent processing elements to the left or top, 
or in the case of processing elements at the top or left boundary of the array, 

20 from operand data formatters. They then separate the Instruction and reformat 
the data to an intemal representation consisting of separate sign, exponent and 
fraction words. This data Is available to the processing element via the X and Y 
intemal data buses. The input registers also compute the sign of the product of 
the two inputs, check for zero operand data and implement the Booth encoder 

25 used during multiplication operations. 

Memory: The memory consists of a Register RIe and a Constant ROM: The 
register file is a 5 word memory used to hold the product (both fraction and 
exponent), accumulator (both fraction and exponent) and temporary results. 
The product and accumulator registers can be swapped under the control of 
30 microcode to facilitate efficient implementation of the pre-alignment operation in 
the floating point addition and accumulation algorithms. The registers can be 
loaded from the R bus, and their contents can be read from either the X or Y 
buses. The Constant ROM stores a number of constants that are used during 
the implementation of the floating point algoritiims. These can be read via the X 
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and Y operand data buses. 

Shifter/Noimalizer: Under microcode control, the shifter can either operate as a 
shifter (for pre-alignment of fractions before addition) or a normalizer. When 
acting as a shifter, it perfonns right-shift operations on one operand datum (The 
5 X operand). The amount by which the datum is shifted is determined by a 

previously computed shift that is applied to the second input to the shifter (the Y 
operand). The shifter can shift 0 to15 bits right within one clock cycle. When 
acting as a normalizer, the Shifter/Normaiizer perfonns either a right shift by 
one bit, or a left shift by 0 to 15 bits within one cycle, in this case the shift is 

1 0 applied to the X operand input to the shifter and is independent of the Y 
operand input. The value of the shift is data dependent. A right-shift is 
performed if the value on the X input is the result of a computation which had 
overflowed (such as in the case of addttion of two normsdized numbers having 
the same exponent). Otiienvise, a left-shift is perfomfied. When acting as a 

1 5 normalizer, the Shifter/Normaiizer at the same time computes the offset 
(exponent offset) ftiat must be applied to the exponent of tiie number being 
normalized in order to compensate for the shift that is applied. Shifting and 
normalization operations that require shifts of greater than 51 bits can be 
implemented by multiple passes through the shifter/nonnalizer. 

20 Arithmetic Unit: The arithmetic unit consists of an input multiplexer, an adder, a 
result shifter and a flags unit. There are two parallel data inputs (X and Y) to the 
arithmetic unit and a single parallel data output (R). The input multiplexer can 
be used to complement and/or left-shift the X operand under control of the 
Booth encoding logic contained in the input registers. This feature is used in the 

25 implementation of multiplication using a modified Booth algorithm. The 
multiplexer can also be controlled directly by the processing element's 
microcode to focilitate the implementation of addition, subtraction and data- 
move operations. 

The adder performs conventional two's complement addition. The carry inpyt to 
30 the adder can be controlled by either the booth encoder logic or the processing 
element's microcode. Both addition and subtraction can be performed by the 
combination of input-multiplexer and adder. Under control of ttie processing 
element's microcode, the output shifter latches either the result of the 
computation or the result divided by 4. This feature is used during partial 
35 multiplication operations. The latched resuft remains valid until the next time the 
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arithmetic unit is used. The latched result can be written onto the resuit bus R. 
A number of flags are set or cleared depending upon the result latohed by the 
arithmetic unlf s output shifter. These flags Include the sign of the result, 
whether or not the result is zero, and whether or not the result Is less than or 
5 equal to 15 (used to support multi-pass shifting during addition pre-alignment). 

Output Registers: The output register module is used to communicate the 
results of computations back from the processing element toward the left 
boundary of an array of processing elements. The ou^ut register can be 
parallel loaded from the arithmetic unit or can be serially loaded from a serial 
1 0 source (often the serial source is from another processing element's output 
register). The output register is unloaded serially. The output register is parallel 
loadable by the arithmetic unit in three parts: sign, exponent and fraction. This 
facilitates conversion from the internal data representation to IEEE 754 floating 
point format 

1 5 During the time when the arithmetic unit is converting the accumulator contents 
into IEEE floating point fomnat. a flag is set to indicate that a register unload is 
in progress (UIP). 

Control and Sequencing: This module includes a microcode ROM, a program 
counter, branch control logic, flags logic, an instruction register, an instruction 

20 decoder, address decoders and timing and control logic. This circuitry is used to 
sequence the processing element through its operations. Each clock cycle, the 
microcode ROM issues a microinstruction to the processing element's datapath 
units, and thereby controls the function and timing of the data operations being 
performed. Data and control flags fed to the branch control togic enable the 

25 processing element to perform data dependent operations required for 
implementation of the floating point algorithms. Relds of the instruction 
transmitted serially to the processing element as part of the {instruction, data} 
2-tuple are also fed to the branch control logic and flags logic of the Control and 
Sequencing Unit These also determine the sequence of microinstructions. 

30 executed by the processing element The instructions specified in the 
{instruction, data} 2-tuple are distinct from the set of microinstructions 
implemented by the processing element The instructions specified in the 
{instruction, data} 2-tuple control the flow of execution of the processing 
element's microcode; 
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The internal data representation used by tlie processing element uses two 32- 
bit data words to represent each IEEE single precision number. One of the two 
words represents the mantissa In 2's complement form, normalized to bit 29. 
The second word represents the exponent using an exponent bias of 229. This 
5 fomiat pro\^des better resolution in the mantissa than IEEE single precision 
format, and the use of a large exponent field virtually guarantees that exponent 
overflow cannot occur. 

Within each processing element, multiplication is feicilitated by the inclusion of a 
modified Booth encoder and multiplexer. The denomnalisation and 
1 0 normalisation operations required by the floating point accumulation or addition 
algorithms are facilitated by ttie repeated application of the shifter circuit which 
can shift up to 15 bits in a single cyde. 

FIG. 2 shows a scalable array processing chip composed of a 5 x 4 rectangular 
array of single precision floating point processing elements which accept serial 
1 5 dataflow operands, and which perform a set of operations on those operands. 
Each operand consists of a 5-bit instruction followed by an IEEE standard 
single precision number. Each processing element is a microcoded ALU with a 
32-bit parallel datapath that includes dedicated hardware ^pport for floating 
point multiplication and addition algorithms. 

20 The array of processing elements is clocked synchronously. The three bit- 
serial links provide communication between processing elements. One link is 
provided for each of the two Input X and Y operands and one for the output, or 
result operand R. As shown in FIG. 2, input data is transfen-ed from left to right 
across the array, and output results are transmitted from right to left Chips can 

25 be cascaded arbitrarily in both X and Y directions. 

The operation of the scalable array processing chip is described with reference 
to the system block diagram shown in FIG. 3. The data interface provides 
communication between the scalable anay chip and the host system. The data 
formatter elements are described separately in a co-pending application 
30 number PL5696 entitled DATA FORMATTER. 

The I/O architecture of each processing element consists of two orthoganal 
data transmission paths for X and Y operands, each consisting of a single one- 
bit delay ceil and a 32-blt data storage register. The X operand path also 
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includes a 5 bit instruction register. Data is input to the array as a sequence of 
{instruction,clata} 2-tupies. These are split into separate Instruction and data 
words on receipt by the input registers. 

Each X data operand consists of a 5-bit Instruction followed by a single 32-bit 
5 IEEE 754 standard floating point number, A variable length gap of several 
clock periods may be present between operands for I/O synchronisation. The 
operand is transmitted in bit serial form into the processing element. When the 
entire {instruction, data} 2-tuple is held within the processing element, it is 
cross-loaded into parallel holding registers. The instruction is decoded and 
1 0 used to control the execution of the floating point algoritiims. The data is 
converted by hardware into the internal extended fomnat. The internal format 
has both extended precision and extended dynamic range when compared with 
the IEEE standard. 

The bit-serial data is bit-skewed on entry to adjacent processing elements on 
15 tiie array boundary. This skew is preserved between adjacent elements witiiin 
the array by passing the data through the single-bit delay stage in each 
processing element before re-transmitting it to the next processing element. 
The use of serial data both minimises the I/O pin count at the array boundary 
and allows adjacent processing elements to both commence and conclude their 
20 computations with a time differential of only one bit period. The advantage of 
the bit-skewing approach over a broadcast architecture is tiiat there is no need 
to drive long buses with large buffers and thereby provides the capability for 
arbitrary expansion of the array. 

Bit skewing has tiie advantage over word-skewing in that less wavefronts are 
25 required to complete a processing task. The bit-skewed approach therefore 
results in the minimisation of job time. The computation time is minimised for 
both a single job and a job stream. 

At ttie completion of a set of computations, an operand wavefront is issued to 
the array which causes the unloading of the results into the output registers of 
30 the processing elements. 

Clocking of the scalable an-ay processing chips is perfomrted by a single phase 
50% duty cycle clock from which all internal timing signals are generated. The 
clock is buffered on entry to the chip and is distributed to each processing 
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element (t is re-buffered within the processing element where it is used as a 
locally synchronous ciocic In addition, each processing element generates a 
second, synchronous clock of the same frequency but with a duty cycle 
detemnined by a self-timed circuit. The secondary clock is used to provide 
5 timing information for bus precharging, data transfers and evaluation of 
execution units. 

FIG. 4 shows schematically the inner-product-step process described by Kung 
and Leiserson. Data is clocked into each processing cell from the left and top 
edges while the results are clocked out from right to left For a matrix product 

10 algorithm, an inner product accumulate algorithm is used in preference to the 
inner product step process common in much of the prior art. The iriner-product- 
accumulate process is depicted schematically in FIG. 5. Data is again clocked 
into the element from the left and top but in this case the result Is formed In 
place. An explicit unload phase is implemented to obtain the result after the 

1 5 computation is complete. An advantage of the inner product accumulate 
algorithm over the inner product step approach is illustrated when matrix 
products are computed for matrix operands which are rectangular. The inner 
product step process requires the recirculation of the result partial product 
matrix. In contrast, the inner product accumulate algorithm computes the result 

20 in-ptace, and Incurs no hardware penalties, irrespective of the length of the 
inner products. 

The sequence of operations performed by the processing elements Is 
determined by the 5 bit instruction transmitted as part of the X operand. The 
five instruction fields and their function are listed in the table below. 

25 



Instruction 


Bit No. 


Function 


ADD 


4 


Floating point add 


LDR 


3 


Convert result to lEE format and load O/P register 


HAD 


2 


Enable result unloading only if active flag set 


SDE 


1 


Set active flag if accumulator contents are non-zero 


CLR 


0 


Clear accumulator prior to computation 



TABLE 1 
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The default operation performed by the PE (i.e., when none of the fields of the 
instruction are asserted) is an inner-product operation implemented as a 
floating point multipiy-accumulate, the input X and Y operands being multiplied 
and accumulated with the contents of the accumulator. 

5 If the CLR field is asserted, the accumulator contents is cleared before the 
computation commences. This generally occurs for the first wavefront of a 
matrix multiplication, and also when executing element-wise operations. The 
accumulator is cleared before the computation is commenced but after the 
ACTIVE flag is set if (the SDE field) is set, and after the accumulator has been 
1 p unloaded into the result register (if the LDR field is set). 

If the SDE field is set AND the value held in the accumulator (from the previous 
operation) is non-zero, an internal flag, ACTIVE (one per processing elemeht) 
is set to indicate that this processing element is an active element Only active 
elements are permitted to unload results during element-wise operations. 

15 If the HAD field is asserted, the operation being performed is deemed to be an 
element-wise (hadamard) operation. If this field is set, only those processing 
elements flagged as active elements (as detemnined by tiieir ACTIVE flags) can 
unload their accumulator contents into their output register R. 

If tile LDR field is asserted, the accumulator contents from tiie previous 
20 computation are converted back to IEEE format and are unloaded into tiie 
processing elements output register. 

During the unloading process, the processing element issues a flag (UIP) to 
indicate that the unload is in progress. 

If the ADD field is asserted, the X and Y operands are added ratiier than 
25 multiplied prior to the result being stored in the accumulator. The HAD and CLR 
fields must also be asserted for matrix addition instructions. 

The element-wise operations of addition and multiplication defined by 

C s A 4* B >A^ere cij ^ aij + bij 

C = A • B where cij = aijbij 
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are performed by first setting the active-element (ACTIVE) flag in a desired 
processing element. This procedure is typically done once during system 
initialization. It is achieved by issuing an instruction with the SDE (set active 
element) field asserted. When this occurs, processing elements that contain 
5 non-zero results in their accumulators set the value of their ACTIVE flag to 
TRUE 

The processing elements accept operand data and return results in IEEE 
standard format. Internally, an extended precision format is used for both the 
mantissa and exponent of the partial results. 

1 0 The internal formats used for the representation of mantissa and exponent are 
as follows: 

2's complement mantissa sgf. fffffffffffffffffffffffffffff 

Exponent Oeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

Exponent Bias 001 000000000000000000000000000000 

15 where: 

s is the sign bit of the mantissa 4-/ is represented as 0/1 
respectively. 

g is a guard bit used to avoid mantissa overflow during 
accumulation. 

20 f is a fraction (mantissa) bit. The mantissa is normalized : the most 

significant fi^action bit is 1 (explicit). 

is the position of the binary point (showing that the mantissa is 

nomnalized). 

e is a bit of the exponent, v\^ich is held in biased form. The 
25 exponent bias is 229. 

If the flags in the anti-active processing elements have been set by a prior SDE 
instruction, and an element-wise multiplication of a matrix A with the unit matrix 
is executed, tiie result of tiie operation is tiie transpose of the matrix A. If an 
arbitrary orthogonal set of elements have their flags set, a pemiutation of the 
30 input matrix will be performed by this element-wise product. 

When an unload (LDR) instruction is received, the accumulator contents are 
converted fiiom the internal format to an IEEE standard form. Numbers outside 
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the range that can be represented by the IEEE single precision format are 
truncated to zero (in the case of results with large negative exponents, including 
IEEE denomnalized numbers) or limited to infinity (in the case of numbers with 
large positive exponents). In both cases, the sign of the zeros or infinities are 
5 retained (unless the result is a true zero, in which case positive zero is always 
returned). 

The IEEE representation of the result is loaded into a separate output register 
which is concatenated with other output registers in adjacent processing 
elements to form an output register chain. The result is output in a serial form 
1 0 through this register chain. 

Matrix algorithms which are elements of the set of primitive operators 
{multiplication, addition, element-wise (or Hadamard) multiplication, 
permutation} are perfomied directly by the processing array. Implementation of 
these operations for operands whose dimension exceeds the size of the array 
15 is possible by mathematically partitioning the operations to a set of operations 
which can be computed separately using the available array size. 

For the particular case where the problem size does not exceed the size of the 
array, recursive algorithms can be implemented which recirculate the output of 
the array bacic to its input^ This can be a useful method to minimise memory 
20 bandwidth requirement in particular applications. 

If a matrix multiplication is commenced with an instruction which does not clear 
the accumulator, the result of the multiplication will be summed with the prior 
result This gives a matrix multiplication/accumulation capability which has 
direct application to the evaluation of complex matrix operations. 

25 FIG. 6 shows the way in which conformal matrix operands are entered into the 
systolic array. Bit-skewing Is indicated by the small offset between adjacent 
rows of A and columns of B. Each element of the processor an-ay computes an 
element cij of ttie result matrix C, by evaluating the inner product 
Ci j = Sk=:0 ^ik^kj • When the last wavefront has been input to the an-ay, 

30 the result matrix may be read from the array. The elements are obtained in tiie 
order shown in FIG. 7. 

If the processing elements on the main diagonal in the array have tineir active 
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element flag set by an arbitrary prior operation as shown in FIG. 8, the 
processor array can be used for the element-wise operations of addition and 
Hadamard muttlpiication. FIG. 8 shows the entry of conformal matrices to a 4 x 
4 subarray of the chip for the purposes of element-wise addition or 
5 multiplication. Only those elements shown as • are used. 

FIG. 9 shows the relationship between rows of data which are output from the 
array after an element*wise operation. Due to the word-length registers present 
in the output register chain, the data is skewed by one word-time plus one bit- 
time. The additional bit-time delay is caused by the bit sicewing of the input 
10 operands. 

In a second embodiment the invention has been implemented In a system 
hosted by a Sun SPARCstation. The matrix processor Is interfaced to the Sun 
SPARGstation via the SBus. This an^angernent is convenient since it ailpws the 
SOAP hardware to operate using virtual addressing, with virtual to physical 

15 translation being perfonmed by the SBus controller in the SPARCstation. The 
host processor and the matrix processor therefore share the same data space, 
so both can interact with the matrix data directly. This approach does however 
have its own disadvantages, the most critical being the fact that the data 
transfer rate across the SBus tends to be quite low due to the overheads of 

20 address translation. 

To compensate for this low data rate, the matrix processor also includes a 
cache memory subsystem. The cache supports burst mode data transfers 
across the SBus on cache misses and can also be used to hold frequently used 
operand matrices (such as coefficient matrices in transform applications) and to 
25 store temporary or Intermediate results. 

A novel cache partitioning scheme has been implemented. The technique 
allows the cache to be dynamically divided into a number of regions that are 
guaranteed not to interact thereby ensuring that fetches for one matrix operand 
do not interfere with fetches for the other. The data controllers determine how 
30 the cache is partitioned on a per-operand/result basis (it is also possible to 
assign a cache partition to the instruction streams) by issuing an 8-bit space 
address along with each address generated. Each bit of the space address can 
be set or cleared, or can take on the value of one of the generated address bits. 
In our system implementation, three bits of this space address are used to 
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control non-cached accesses, temporary matrix accesses and temporary matrix 
initialization. Four bits are used to partition the cache into up-to 16 independent 
regions. Use of the temporary matrix control bits of the space address allows 
temporary result matrices to be stored entirely within the cache without being 
5 written out to the host. In fact, such matrices are entirely in>^slble to the host 
processor. The maximum data throughput obtainable using the cache is 12.5 
Mwords/second. 

The two custom chips implemented during the development of this system are 
a processing element array chip and a data controller chip. Both chips were 
1 0 designed using a generic 1 .2 micron double layer metal CMOS process rule-set 
and were retargetted for fabrication using a 1.0 micron process using a gate 
shrink. 

The processing element an^ay chips are full custom integrated circuits each 
containing an array of 4 rows by 5 columns of floating point processing 
1 5 elements. Because the overall computation rate Is limited by the available data 
bandwidth, tiie speed of computation of the processing elements If not overiy 
important. Therefore, the architecture has been designed to yield processing 
elements (PEs) that are physically small rather than being particulariy fast. 
Each complete floating point unit occupies only 2.7sq mm. 

20 The processing element does not include a dedicated hardware multiplier, but 
is implemented as a simple microprogrammed 32-bit datapatii with hardware 
support to aid the floating point computations, as illustiBted in FIG. 5. 

The PE hardware incorporates a booth encoder and multiplexer to facilitate 
multiplication using an iterative modified bootii algorithm, and also a 

25 shifter/hormalizer that can be used for pre-addltion alignment as well as post 
addition nomialization. When used as a normalizer, the shifter has the ability to 
compute the amount by which the exponent must be adjusted during the same 
time that the nomialization occurs. Computation of the floating point arithmetic 
operations (multiply/accumulate, multiply or add/subtract) are completed within 

30 40 clock cycles. 

The processing element array chip accepts IEEE single precision floating point 
numbers as inputs and feeds results back through the data controllers in the 
same format Intemally, a proprietary number representation is used, including 
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a 31 bit exponent that virtually eliminates the possibility of exponent overflow. 

The chips operate at 20MHz clock speed, achieving around 20 MFLOPS peak 
performance per chip. Processing anays of arbitrary size can be buiit with no 
external components simply by stacking the chips to fomfi a two dimensional 
5 array. The pin*out of the chip is such that 1-to-1 connection of inputs and 
outputs of adjacent chips can be made. All communication to and from the 
array is via the edge elements of the array. Operand data enters the array on 
left and top edges. This data is known as the X and Y operand data 
respectively. The result data (R) emerges from the left edge of the array and 
10 can be extracted independently from the application of operand wavefronts 
(that is, the operand and result streams operate in parallel). 

The only global signals in the array are clock and reset Because ail 
communication is local (nearest neighbour only), the system Is insensitive to 
clock skew from one side of the array to the other. The only requirement is that 
15 the skew between adjacent PEs is kept under control. This can be readily 
achieved by orderly layout of clock routing and/or insertion of clock buffering. 

The processing elements are low power devices due to their architecture. The 
entire chip containing 20 processing elements dissipates less than half a watt. 
This corresponds to less than 5mA per processing element at 20MH2 
20 operation, or 5mA per MFLOP. 



Number of Transistors 


270 000 


Die Size (Pad to Pad) 


8.56mm x 8.35mm 


Transistor Density 


3800 T/sq mm 


Power Dissipation 


0.5 Watts 


Package 


68 CLCC 


Floating-point Performance 


20 MFLOPS @ 20 MHz 


Design Style 


Full Custom 



TABLE 2 



The performance attained by the apparatus of the second embodiment for a 
range of applications is shown in Table 3. 
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Applicatton 


Exec. Time 


Performance 


3450 point 1 D Fourier transform 
using 2D factorization 


20 msec 


130 MFLOPS 


2D Fourier Transfonn of 380380 
point image 


1385 msec 


66 MFLOPS 


4000 tap FIR RIter 


35 msec per 1000 
data samples 


210 MFLOPS 


10th order Matrix polynomial 
evaluation of 60 x 60 complex 
matrix 


136 msec 


114 MFLOPS 


QR factorization of 59 x 60 Matrix 


561 msec 


87 MFLOPS 



TABLE 3 
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1 . A processing element suitable for use in a scalable array processor 
comprising of at least one input register means adapted to receive and process 
serial operands in the form of {instruction, data} 2*tupies, a memory means 
adapted to store temporary results and constants, a computing means adapted 
5 to perform iogical operations, an output register means adapted to output 
results from the processing element, a control and sequencing means adapted 
to control the operation of the processing element; 
a plurality of data buses adapted to provide communication between the 
plurality of means. 

10 2. An apparatus as in claim 1 wherein the computing logical means 

consists of a shifter/normaliser means adapted to shift/normalise data and an 
arithmetic means adapted to perfoni) logical operations such as but not limited 
to addition, subtraction and partial multiplication operations. 

3. An apparatus as in claim 1 wherein the processing element is adapted 
15 to perform floating point multiply, floating point add and floating point multlply- 

accumulate which can be used for inner product accumulate operations. 

4. An apparatus as in claim 1 wherein the input register Is adapted to output 
a copy of an input operand bit with a one clock period delay. 

5. An apparatus as in claim 1 wherein there are N input registers and the 
20 processing element Is suitable for use in a N-dimensional scalable array 

processor. 

6. An apparatus as in claim 5 wherein N can be any positive integer. 

7. An apparatus as in claim 1 wherein the input registers are adapted to 
convert the input serial data to an internal representation comprising separate 

25 sign, fraction and exponent. 

8. An apparatus as in claim 1 wherein the memory means consists of read 
only memory for storage of constants and a read/write memory for storage of 
temporary results. 

9. An apparatus as in claim 2 wherein the shifter/normalizer means is 
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adapted to perform binary weighted barrel sliifting wherein the shifter function is 
determined by a control input to the shifter/normalizer and the normallzer 
function effects a data dependent shift of up to 15 bits within a single cloclc 
cycle. 

5 10. An apparatus as in claim 2 wherein the arithmetic means implements 
logical operations such as but not limited to floating point addition, multiplication 
and multiply-accumulate algorithms using a parallel microcoded data path. 

11. An apparatus as in claim 2 wherein the arithmetic means comprises a 
logical unit such as but not limited to an input-multiplexer, an adder, an output 

1 0 shifter, a flags unit and a control unit. 

12. An apparatus as in claim 1 wherein the output register is adapted to be 
loaded in parts to enable the conversion from the intemal representation to 
IEEE 754 floating point format and the output register is adapted to be parallel 
loaded from the arithmetic means or serially loaded from a serial source and 

15 the register is unloaded serially. 

13. An apparatus as in claim 1 wtierein the control and sequencing means 
includes timing and control logic, a microcode ROM, address decoders, branch 
control logic, flags logic, instruction register, instaiction decoder and a program 
counter. 

20 14. An apparatus as in claim 1 wherein there are three data buses, an X bus, 
a Y bus and a R bus and where the X and Y buses are called operand buses 
and the R bus is called the result bus. 

15. An apparatus as in claim 1 wherein there is an accumulator comparison 
means. 

25 16. A scalable array processor chip comprising an array of processing 
elements each said element including at least one input register means 
adapted to receive and process serial operands in the form of {instnjction, data} 
2-tuples, a memory means adapted to store temporary results and constants, a 
shifter/normalizer means adapted to shift or nomnalize data, an arithmetic 

30 means adapted to perform logical operations such as but not limited to addition, 
subtraction and partial multiplication operations, an output register means 
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adapted to output results from the processing element, a control and 
sequencing means adapted to control the operation of the processing element, 
a plurality of data buses adapted to provide communication between the 
plurality of means and wherein each processing element has means for 
5 communication only vvnth adjacent elements. 

17. An apparatus as in claim 16 wherein the array of elements comprise an 
interconnected lattice of at least one dimension. 

18. An apparatus as in claim 16 wherein the array of elements comprise an 
interconnected lattice of at least two dimensions. 

10 19. An apparatus as in claim 1 6 wherein the scalable array processor chip is 
adapted to perform at least the functions of computing the product of two or 
more matrices, computing the element-wise product of two or more matrices, 
computing the sum of two or more matrices, penmuting the rows and columns of 
a matrix and transposing a matrix. 

15 20. A computing apparatus comprising a host processor, at least one 
scalable array processor chip and a plurality of data formatters wherein the 
scalable array processor chip(s) and plurality of data formatters are adapted to 
perform matrix operations otherwise performed by the host processor. 

21 . An apparatus as in claim 20 wherein the apparatus includes a memory 
20 cache adapted to store operand data and temporary or intermediate results. 

22. A method of performing matrix operations comprising the steps of : 

(a) providing a plurality of processing elements in the form of an an^ay 
adapted to perform systolic processing operations; 

(b) receiving operand matrix data for processing from a host or data source; 
25 (c) formatting the operand matrix data in a data formatter by adding an 

instruction to form an {instruction, data) 2-tuple; 

(d) transferring sets of 2-tuples to the processing element array to cause the 
processing elements to process the data in accordance with the instruction; 

(e) repeating the steps (b) to (d) a number of times said number of times 
30 being dictated by the matrix operation being performed; 

(f ) unloading the results of tiie matrix operations into an output result 
register of the processing elements (under control of an instnjction specified by 
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the operand 2-tuple); 

(g) transferring the contents of the output result registers held within the 
plurality of processing elements back to data formatters as result wavefronts; 

(h) storing the result wavefront data back to a host or data sink; and 

5 (i) repeating the steps (f) to (h) a number of times said number of times 
being dictated by the matrix operation being performed. 
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